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ABSTRACT 


Geometric and variational techniques, along with the method of Lagrange 
multipliers, optimal control theory, and elementary calculus are the tools used to 
derive some generalizations of Snell's Law. The focus is both on Snell’s Law of 
Refraction and, on the reflection principle under various conditions. 
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I. INTRODUCTION 


Snell’s Law of refraction and the principle of reflection have both been known 
for several hundred years. The intent of this paper is to develop some 
generalizations of these observed phenomena. We derive these generalizations by 
imposing certain constraints upon the problem and show that in all cases our general 
results reproduce the now familiar forms of the Law of Refraction and of the 
reflection principle when suitably restricted. 

In this paper we make use of variational methods along with techniques from 
optimal control theory to develop a generalization of Snell’s Law of Refraction in 
which the velocity of light in the two media depends upon direction alone. We also 
develop several interesting inverse problems which involve the determination of 
various indices of refraction corresponding to certain families of curves. For 
instance, if we desire to move from a point A to a point B along a cycloid, in order 
to avoid an obstruction in our path, we can determine the necessary index of 
refraction needed to do so. 

"The Dutch astronomer and mathematician Willebrord Snell (1591-1626) 
discovered the law of refraction by considering a ray of light that passes from one 
medium, in which light travels at a velocity v,, to a second medium, in which light 


travels at a velocity v,." [Ref. 1]. 








Fermat’s Principle states that the time elapsed in the passage of light between 
any two fixed points of a medium is an extremum with respect to possible paths 
connecting the points [Ref. 2]. Since the distance travelled is proportional to the 
time within each medium, the path of minimum time in a medium which has a 
constant velocity of propagation is a straight line. If however, we are given two 
media in contact with propagation velocities v, and v, (v, ~ v,), the path followed 
between a point (X,,Y,) in medium 1 and a point (X,,Y,) in medium 2 consists of two 
line segments meeting at a point (X,,Y,) on the interface of the two media. This 
"bending" of the light ray is called refraction, and hence, Snell's Law of Refraction, 
which is given by: 


sin(,) _sin($,) 


vy Vo 


where ®, and ®, are the angles formed by the ray and the normal to the interface. 
The paths in which light particles travel are represented by rays, which indicate the 
direction of propagation. These rays are perpendicular to the wave fronts, regardless 
of isotropy. An isotropic medium is one in which the properties of the medium are 
the same in all directions from any given point. A medium is homogeneous if its 
properties do not vary from point to point. Wavefronts can have many shapes. If 
the disturbances are propagated in a single direction, the waves are called plane 
waves. "At a given instant, conditions are the same everywhere on any plane 


perpendicular to the direction of propagation.” [Ref. 1:p. 301). 








It has been observed that if a plane light wave falls on another plane surface 
the light beam is reflected from the surface (Ref. 1]. The angle of incidence, ®,, and 
the angle of reflection, ®,, measured from the normal to the surface satisfy the 
principle of reflection which is given by: 

®, = 9, 

In the study of acoustic wave propagation, Snell's Law of Refraction and the 
law governing the reflection principle are usvally derived from the reduced 
Helmholtz equation. Although this thesis concerns generalizations of these laws 
using variational methods, a brief outline of this process is included to illustrate 
other available techniques. 

Figure 1.1 depicts plane wave reflection and refraction from a planar interface. 
The equation for the incident plane wave, P, is given by: 

P= eked, __, k,(xsin6,- yoas®,) 


where the incident amplitude has been normalized to 1, the vector ¥ is a position 


vector, and d is the vector which points in the direction of propagation of the wave. 


Using subscripts u for the upper medium, | for the lower medium, i for incident, r 
for reflected, and t for transmitted we have ‘he folicwing system of equations: 
VP, +k’,P, = 0, y > 0 


VP, +k’,P, = Oy > 0 


VP, +k*P, = 0, y < 0 









UPPER MEDIUM 


LOWER MEDIUM 


Figure 1.1 Plane Wave Reflection and Refraction 


-iwt 


where k =w/c = wave number. Ane“ time dependence is assumed throughout the 
following analysis. Allowing for continuity of pressure at the interface, we have the 
condition: 


P,=P,+P,=P,=P, 


or, 








P+ PP, 
For an ideal fluid, we have: 


VP=-pVob=iwpi 


where w is the fluid velocity. After taking the dot product of both sides of this 
equation with # we have: 
h-VP=i -iwpit 


where # is the unit normal to the interface. This equation yields: 


OP. 
— =iwpv, 
where v is the fluid velocity normal to the interface. Our second boundary condition 


at the interface equates normal! fluid velocities, yielding: 
1a 1 oO 
—£(P,+P,)-—<(P), 
Py oy : ? P; 3, : 


where ¢, and ,¢, are densities of the two fluid mediums. We define reflected and 
transmitted plane waves as 


P =Re**4 


and 


P,=Te aed 


Substituting these into the two interface boundary conditions produces: 








ck sin8; ik xsin8, ik, xsin6, 
etkemn8, | ppikersin®, _ or, ik,ssin8, 


k yee -k ror 
—*|-cos0,e""="° +cos0,Re"™"*'|- —“cos6,Te #9 ie 
Py Py 

In order for these two equations to be satisfied for all x, the phases must be equal 
along the interface. From this we have: 


ik,sin©O, = ik,sinO, = iksin®,. 


This equation yields both the law of refraction and the law governing reflection: 





sin8, sin0, 

C, C, 
and, 

0, = O.. 


The solution to the problem for R and T is that 


cos8, - a Bcos8, 
 cos6,+a Bcos8, 


2cos8, P, C, 
=—______—_—.,, where a = — and B=—. 
cos0; +a Bcos6, P, C, 


Special cases such as Brewster’s angle and the critical angle where total internal 
reflection occurs can be studied in the context of how R and T vary for differing 
material properties and indices of reflection. 

Reflection and refraction are both important topics in the study of acoustic, 
electromagnetic, and elastic wave propagation. The following list of books all have 
good treatments of this subject matter: Halliday and Resnick [Ref. 1], Jones [Ref. 


3}, Lorrain and Carson [Ref. 4], and Achenbach [Ref. 5]. 





Il. MATHEMATICAL PRELIMINARIES 


A. INTRODUCTION 

The purpose of this chapter is to present a brief overview of the basic 
mathematical concepts which are employed in future chapters. The important 
concepts include the idea of a functional, the simplest variational problem, the 
Euler-Lagrange necessary condition for the solution of the simplest variational 
problem, special cases of the Euler-Lagrange equation, an introduction to optimal 
control, Pontryagin’s Maximum Principle, and, the use of Lagrange multipliers. The 
following references were used in the development of this material: Gelfand and 
Fomin [Ref. 6], Weinstock [Ref. 2], Leitmann [Ref. 7}, Ewing [Ref. 8], Weir [Ref. 9], 
Bliss [Ref. 10], Pontryagin [Ref. 11], Taylor and Mann [Ref. 12], and, Fleming and 


Rishel (Ref. 13]. 


B. FUNCTIONALS 

Variable quantities called functionals play an important role in several areas 
of mathematics, physics and mechanics. A functional is an operator which maps a 
function onto a real number. 

For the problem of wave propagation in a host medium, consider all possible 
paths joining two given points A and B in the plane. “Suppose that a particle can 


move along any of these paths, and let the particle have velocity v(x,y) at the point 





(x,y). We obtain a functional by associating with each path the time required for the 


particle to traverse the path." [Ref. 6:p. 1] Theoretical proofs of some results in 
calculus of variations depend upon concepts from functional analysis. For 
applications of functional analysis and variational methods in engineering, see Reddy 
{Ref. 14]. For a discussion on applications of functional analysis to time optimal 


control problems, consult Hermes and LaSalle [Ref. 15]. 


C. THE SIMPLEST VARIATIONAL PROBLEM 
The simplest variational problem invclves the determination of the maxima or 
minima of a functional of the form: 
So=f"Fesy,y dx 
with fixed boundary conditions: y(a) = A and y(b) = B. An example of the simplest 
problem is finding the shortest plane curve joining two points A and B. That is, the 


problem is to find the curve, y = y(x). for which the functional 


b 2 
Iyl=f Vi+y"ax 
a 
achieves its minimum. The curve in question turns out to be a straight line 


connecting the two points A and B. 


D. THE EULER-LAGRANGE NECESSARY CONDITION 
If there exists a continuous, twice-differentiable function that minimizes the 
simplest problem, what then is the differential equation which satisfies this function? 


The answer to this question leads to the implementation of the Euler-Lagrange 








equation. We include the following sketch of a derivation of the Euler-Lagrange 


equation fer completeness. For an in-depth development of this matter, see 
Weinstock [Ref. 2:pp. 20-22]. We know that f(x+h) - f(x) = f’(x)h + higher order 
terms in h. We apply this to our simplest variational problem by considering two 
functions y and y+h satisfying the same boundary conditions. Let y(a)= A and 
(y+h)(a) = A, then h(a) =0 The same argument holds true for y(b) = B. Now, 
define the following: 


AJ = Jfy+h} - Jy] 


AJ= IF o.y+hy’+h')-Foay,y Vax 


and use a Taylor series expansion to arrive at: 


AJ= if °(F h+F,h))dx + higher order terms. We also define dJ[h] to be the first 
a 


variation: 


SJ[h] = ibe (F, h+F,h))de. 


Setting the first variation equal to zero, and integrating by parts we obtain: 
b b b b fb, d 
[ Fyhdes [OF hides | UF ade Fh \¢- [he (Fide = 0. 
Now, since h(a) = h(b) = 0, 


b 
Fish |,=0 . 


This leaves us with: 








ffle,-ZF | hdx = 0. 
a dx y'| 
By the fundamental lemma of the calculus of variations, (see Gelfand and Fomin 


(Ref. 6:p. 9]), it follows that: 


if, f "a(x)h(x)dx=0, where h(a) = h(b) = 0, and a and h are C[a,b] then a(x) = 0 


for all x € [a,b]. Then, 
d 
ee =Q. 
Thus, the Euler-Lagrange equation for the simplest problem is: 
F,-2(,) =0. 
Thus, a function y{x) which satisfies the Euler-Lagrange equation and the associated 
boundary conditions is by definition the extremizing function. 

The use of the Euler-Lagrange equation does not imply that the solution is an 
extremuia. That is. a solution to the Euler equation presupposes the existence of a 
maxima or a minima. Just as in the ordinary differential calculus, if we consider 
f(x) = x’, setting f(x) = 3x? = 0 yields neither a minima or a maxima. Another 
necessary condition is the Legendre necessary condition which states that for the 


functional: 


10 





Jty]= i, FQyy dx, 


y(a) = A, y(b) = B to have a minimum for the curve y = y(x), the inequality 
F,, 2 0 must be satisfied at every point of the curve. Other necessary conditions 


are further developed in Leitmann [Ref. 7], and, Bliss [Ref. 10}. 


E. SPECIAL CASES OF THE EULER-LAGRANGE EQUATION 

The Euler-Lagrange equation plays an important role in the calculus of 
variations and is normally a second-order differential equation. There are some 
special cases for which Euler’s equation can be reduced to a first-order differential 
equation. For instance, if the functional has no y dependence Euler’s equation 
becomes F,. = C. If the functional has no x dependence then the corresponding 
Euler equation is F - y’F,. = C. Finally, if the functional does not depend upon y’, 
Euler’s equation reduces to F,(x,y) = 0 and therefore is not a differential equation. 
For further treatment of this subject, see Gelfand and Fomin [Ref. 6:pp. 18-19], and, 


Weinstock [Ref. 2:p. $2]. 


F. OPTIMAL CONTROL 

The problem of optimal control is intimately related to certain problems in the 
calculus of variations. The following is an outline of this process and is taken from 
Gelfand and Fomin [Ref. 6:pp. 218-219]. Suppose the state of a physical system is 


characterized by n real numbers x’,...,x*, forming a vector x = (x’...,x") in the 














n-dimensional phase space X of the system. Further, suppose the state varies with 


time in the way described by the system of differential equations: 


i + 
Sf ok "uu \(i=1,....). Here the k real numbers u’,...,u" form a vector 


u = (u’,...,u') belonging to some fixed control region 22, which we take to be a subset 
of k-dimensional Euclidean space. The f(x,u) are n continuous functions defined for 
ail x « X and all u € Q. Now suppose we specify a vector function u(t), t, ¢ t S t, 
called the control function. with values in Q. Then, substituting u = u(t) in the 


differential equations above yields: 


a =f if} x" l(),...4 (ONG = 1,...9). 
For every initial value of x = x(t,), this system has a definite solution, called a 
trajectory. The aggregate, U = {u(t),tot,Xo}, consisting of a control function u(t), an 
interval (t,t;] and an initial value x, = x(t), will be called a contro] process. Thus. 


to every control process, there corresponds a trajectory. Now let f°(x'.....x"u'....,U") 





be a function which is defined, together with its partial derivatives, of (i=1,....%), 
ox! 


for all x « X, and all u € Q. 


To every control process U, we assign the number ip " £(x,wdt. Then, the 
t 


control process is said to be optimal if the inequality. J[U] < J[U*] holds for any 








other control process U* carrying the given point x, into the point x,. By the optimal 
trajectory, we mean the one corresponding to the optimal control process. The 
following books have detailed explanations of this material and are recommended 
for further study: Pontryagin {Ref. 11], Leitmann (Ref. 7:pp. 79-83], and Berkovitz 


[Ref. 16}. 


G. PONTRYAGIN’S MAXIMUM PRINCIPLE 

We now state the main theorem of optimal control theory, the maximum 
principle. For further reading on this theorem and its proof, see Pontryagin [Ref. 
li:p. 11], Gelfand and Fomin [Ref. 6:p. 222], and Leitmann [Ref. 7:p. 118]. 

Let U = {u(t),tyt,.x 9} be an admissible control process, and let x(t) be the 


corresponding integral curve of the system, 


Sfx = (),1,....n), passing through the point (0.x,',....%.") for t = 0, and 


satisfying the conditions: x'(t,) = x,’.....x°(t,) = x," fort = t, Then if the control 
process U is optimal there exists a continuous vector function W(t) = (%(t). 


¥,(t).....8,(t)) such that: 


avy, = 
1. The function v(t) satisfies the system, ——=-)}> PCH y 


G = O.1...n). for x = x(t) u = u(t). 

















n 
2. For all t in [t,t,], the function H(¥,x,u)= > P f*(x,u) achieves its maximum 
a=0 
for u = u(t). 
3. The relations, ¥%(t,) < 0, G[¥(t,),u(t,)] = 0 hold at the time t, where 


Gy) = SYP HR x1). 


H. LAGRANGE MULTIPLIERS 
"A necessary condition for a minimum (or maximum) of F(x,y,....z) with respect 
to variables x,y....,z that satisfy G(x,y,....z) = C, (i = 1,2,....N), where C, are given 


constants is: 


“ox ay 
N 
where F'= F+ > 4,G, The constants ,, Ax... Ax. introduced as undetermined 
C=1 


Lagrange multipliers, are evaluated together with the minimizing (or maximizing) 
values of x,y,...,2 by means of the set of equations consisting of those listed above.” 
[Ref. 2:p. 6]. 

Taylor and Mann [Ref. 12], and Weir (Ref. 9] give good treatments of the 


Lagrange multiplier method with several examples. 








Il. REFLECTION 


A. INTRODUCTION 

Mathematical problems are frequently inspired through observations of various 
physical phenomena. In fact several solutions to mathematical problems are inspired 
by nature. The physics of the problem provide us the clues, without which we would 
have very little chance of succeeding in finding solutions. One of those problems 
suggested by nature involves the reflection of light particles on a plane mirror. 

We know from Fermat’s Principle that a light particle will choose the path 
which will get it to its destination in the shortest amount of time. We also know that 
the shortest distance between two points is a straight line. But what happens to 
these light particles when they undergo a reflection? Do they still choose the 
shortest path and if so, what is the shortest path? This chapter answers these 
questions and others. The material in this chapter, although well known, is 
presented for completeness. We review the classical results of reflection and 
complete our discussion with a derivation of the reflection principle using the 
method of Lagrange multipliers. The following references were used extensively in 
the preparation of this chapter: Courant and Robbins [Ref. 17], Polya [Ref. 18], 


Kazarinoff (Ref. 19], Williams [Ref. 20], and Chakerian and Ghandehari [Ref. 21]. 












B. HERON’S PRINCIPLE 
Our problem, as shown in Figure 3.1 is as follows: given two points A and B, 
both on the same side of a straight line L, all in the same plane, what is the point 


C on L which minimizes the sum of the distances from A and B? 





Figure 3.1 Heron’s Problem 


That is, find the point C such that: AC + BC is a minimum. If we can show that: 


AC’ + BC’ > AC + BC 





then we prove that the point C does in faci minimize the sum of the distances from 
A and Bto C. We know that: 
AC’ + BC’ = AC’ + B’C’ 
for B’ the point reflected in the line L through the perpendicular to L at E. This 
follows since the triangles BEC’ and B’EC’ are congruent triangles and the line L 
is the perpendicular bisector of the line segment BB’. Now by the triangle 
inequalitv, we have: 
AC’ + B’C’ > AB’ 
Looking at Figure 3.1, we can see that: 
AB’ = AC + CB’ 
Now by the same argument used above (congruent triang!es), we know that: 
CB’ = CB. 
Then it follows that : 
AB’ = AC + CB. 
Therefore, 
AC’ + BC’ > AC + BC 
as we desired. So we have in fact established that the point C does minimize the 
sum of the distances from A and B. It follows from geometry that the angles ®, and 
@, are equal and thus we have the law which governs the principle of reflection: 
®, = ®, 
Therefore, the reflected ray of light takes the shortest possible course between 


the two points A and B. This reflected ray stays within the plane of incidence since 
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any other possible direction of reflection would not produce the shortest possible 


path from the point A to the point B. This discovery is due to Heron, the 
Alexandrian scientist of the first century A.D. [Ref. 17]. "Simple as it is, Heron's 
discovery deserves a place in the history of science. It is the first example of the use 
of a minimum principle describing a physical phenomenon.” [Ref. 18] 

Heron’s reflection principle can be generalized in Minkowski spaces. A 
Minkowski space is simply a finite dimensional normed linear space. A 
generalization of Heron’s reflection principle using Lagrange multipliers in a 
Minkowski plane and related applications to conics is given in Ghandehari [Ref. 22]. 
For a generalization of Steiner’s problem (Fermat’s problem) see Chakerian and 
Ghandehari [Ref. 21]. 

Consider the following problem: given an area A and one side of a triangle. c, 
which triangle is the one for which the sum of the other sides a and b is the 
smallest? See Figure 3.2. [Ref. 17:p. 332] 

Prescribing the side c and the area A of a triangle is equivalent to prescribing 
the side c and the height h onc, since A = ¥(ch). Therefore, our problem is 
reduced to finding the point R such that the distance from R to the line PQ is equal 
to the height, h and such that the sum a + bis a minimum. From the first condition, 
it follows that R must lie on a line parallel to the line PQ at a distance h. This R 
is given by Heron's Theorem for the special case where P and Q are equally distant 


from L. 
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Figure 3.2 Application of Heron’s Principle 


C. STEINER’S PROBLEM 

The following example comes from Polya [Ref. 18:p. 145], and illustrates the 
use of the reflection principle. It is commonly known as Steiner’s problem named 
after Jacob Steiner, the famous representative of geometry at the University of 
Berlin in the early nineteenth century. Three towns intend to construct three roads 


to a common traffic center which should be chosen so that the total cost of the road 
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construction is minimized, assuming cost is proportional to length alone. Simply put, 
given three points, find a fourth such that the sum of the distances from the three 
points to the fourth is minimized. See Figure 3.3. Let A, B, and C denote the three 
towns and X the fourth point in question. Also, assume that all interior angles a, b, 


and c are all Jess than 120°. 


Figure 3.3 Steiner’s Problem 
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We want to minimize the following expression: C = AX + BX + CX. If we fix one 


of these distances, say CX, then we have a problem which resembles Heron’s 
problem. Therefore, we let CX = r and minimize AX + BX. In this case, we are 
dealing with reflection in a circular mirror. We know that light particles move in 
such a way that the angle of incidence is equal to the angle of reflection. So in our 
problem the angle <AXB must be bisected by a straight line passing through C and 
X. Similarly the angles <AXC and <BCX must also be bisected by lines through 
B and A respectively. As such, the three straight lines joining X to A , B, and C 
dissect the plane into six angles, the common vertex being X. It then follows that 
all six angles are equal and equal to 60. Therefore the three roads diverging from 
the traffic center do so at 120° angles. If one of the interior angles is at least 120°, 
then the common vertex which minimizes the sum of the distances from the three 


points A, B, C will be that vertex which has an interior angle of at least 120°. 


D. REFLECTION PROPERTIES OF CONICS 


1. Ellipse 
The reflection principle applied to an ellipse indicates that the focal radii 
make equal angles with the tangent line. The following derivation is taken from 
Kazarinoff [Ref. 19]. Let L be a tangent to an ellipse with foci F, and F,, see Figure 
3.4. Let P be the point of tangency, and let Q be any other point on L. Since Q is 


outside of the ellipse: 
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QF, + QF, > PF, + PF,. 
Therefore, F,PF, is the shortest path from F, to L to F, and the angles a and B are 


the same. 





Figure 3.4 Reflection Properties of an Ellipse 











2. Parabola 
An important property of parabolas is that any ray of light emanating from 
the focus of a parabolic mirror is always reflected along a line parallel to the axis 
of the parabola. This property is equivalent to saying that the angles a and B, see 
Figure 3.5, formed by lines through the focus F and an arbitrary point P on the 
parabola and through P parallel to the axis of the parabola are equal. For a proof 


of this property, refer to Williams [Ref. 20]. 





Figure 3.5 Reflection Properties of a Parabola 
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3. Hyperbola 
The tangent properties of the ellipse are related to those of the hyperbola. 
Figure 3.6 depicts a hyperbola with P and Q as its foci. The tangent line L passing 
through an arbitrary point R on the hyperbola bisects the angle subtended at that 
point by the foci of the hyperbola. [Ref. 17:p. 335] Therefore, the angles a and 8 


are the same. 





Figure 3.6 Reflection Properties of a Hyperbola 





E. A GENERALIZATION OF REFLECTION 

We now turn our attention to a derivation of the law governing the principle 
of reflection involving a ray of light passing through a homogeneous medium, see 
Figure 3.7. We make use of the method of Lagrange multipliers in developing our 


solution. 





Figure 3.7 Reflection in a Homogeneous Isotropic Medium 





We know that: 





Similarly, for ®,: 
OB = bsec(®,) 
Applying Fermat's Principle. we wish to minimize the following expression: 
asec(®,) + bsec(®,), subject to the constraint that: 
atan(®,) + btan(®,) = c. 
Now construct a new function, call it h such that: 
h( 5,4) =asec(o,) +bsec(,)-Afatan(,)+btan(,) -c] 
Applying the method of Lagrange multipliers [Ref. 9], we set: 


which yields the following system of equations: 
oh 2 
anes aed a $,=0 
1 


Dh geect,tand, - Absec*, = 0 
0; 


which yields: 
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asing, = Aa 





cos’, cos’, 
bind, Ab 
cos), cos’, 


-atang, - Hand, -c =0 


We then have: 


sing, =A and sing,=A 


which yields the law governing the principle of reflection: 


9, > 4. 











IV. REFRACTION 


A. INTRODUCTION 

"A stick bobbing in the water looks sharply bent. We reason that the light that 
follows a straight course in the water as in the air undergoes an abrupt change of 
direction in emerging from the water into the air." (Ref. 18] | This is the 
phenomenon of refraction. The light that travels with a known velocity through air, 
travels with a different velocity through water. Such a difference in velocity explains 
in part the phenomenon of refraction. This chapter will deal explicitly with Snell's 
Law of Refraction and is the focus of this paper. We develop several derivations, 
the most important of which is one in which the velocity of light depends upon the 
direction the light takes in each medium. We use the method of Lagrange 
multipliers and optimal control theory to derive this generalization. The following 
references were useful in the preparation of this chapter, Apostol [Ref. 23]. Polya 
[Ref. 18], Weinstock (Ref. 2], Taylor and Mann [Ref. 12], Pontryagin [Ref. 11], 


Gelfand and Fomin [Ref. 6], and, Halliday and Resnick [Ref. 1]. 
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B. A SIMPLE EXAMPLE 

Our first development, as an introduction, comes from A. V. Baez, which was 
published by Apostol [Ref. 23]. 

It is known that of all possible paths connecting two points A and B in space, 
a light ray leaving A chooses the path which enables it to reach B in the least time. 
This is Fermat's Principle from which the law of refraction and the reflection 
principle follow. 

Consider the following problem: A man is in a boat at a point P, one mile from 
the nearest point A on the shore. He wishes to go to a point B which is farther 
down the shore, M miles from A. He also knows that he can walk faster than he 
can row the boat. If he can row r miles an hour and walk w miles an hour, toward 
what point C should he row in order to reach B in the least amount of time? See 
Figure 4.1. Looking at Figure 4.1 if we let AC = x, Snell’s Law of Refraction 


indicates that: 


sind, + 


sind, w 


Let ®, = 90° and use geometry to obtain: 





Solving for x we have our solution: 





1 
i) 
i) 
{ 
! 
1 
i] 
1 
4 
4 
4 
i} 
{ 
i 
f) 
t 
' 
| 
1 
' 
1 





WATER 


M Miles ———-4B LAND 





Figure 4.1 Application of Law of Refraction 


This example serves to illustrate the fact that we can analyze and solve certain 


problems avoiding the use of calculus by simply applying Snell’s Law of Refraction. 
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C. GENERALIZATIONS OF SNELL’S LAW OF REFRACTION 


1. Optically Isotropic Media 
Our first derivation of Snell’s Law of Refraction involves the passage of 
light between two points in two optically isotropic media. See Figure 4.2. That is, 
the velocities of the particles in both media are constant but different [Ref. 2]. 
Consider for example the path of a light particle emanating from a source at (X,,Y,) 


in water through a thin plate of glass, to a point (X,,Y,) also in water. 


(X1,¥1) 








——_ Incident 


Reflected ~~ Ray 


Ray 


Refracted 


Ray XN 


Figure 4.2 Refracted, Reflected, and Incident Rays 








Applying Fermat’s Principle from the point (X,,Y,) in medium | to a point 


(X,,Y;) in medium 2, with respective light velocities of v, and v, we arrive at: 


T V%, -X,) +, -Y,y “ eS -x,y a ie yy 


vy vz 


where T = time = arc length/velocity = > <-> dt 


dt 


Again Fermat’s Principle claims that the light particle will take the path of shortest 
time between two points or, 
aT X, -X, 7 X,-Xp -0 
aX, Voy (Xp -Xp)° +(¥q- Yo)? vyy(X, - XQ)? + (Y, ~ Yo)? 
X, -X, X,- 
1 0 : and sind, = 0 X, 
J (X, -Xo)?+(¥, -yp)" y(X-X,) +(¥o-¥,)° 


Then the expression becomes: 


Since sing, = 


as 


vy Vv, 


2. Optically Isotropic Media Using Lagrange Multipliers 
We now attempt to solve the same problem, that is the propagation of 
light in two isotropic media using the method of Lagrange multipliers [Ref. 12]. See 


Figure 4.3. 


Looking at Figure 4.3 it is easy to see that cos(@,) 7 or, 
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OA =asec (®,). Now applying Fermat’s Principle we want to sum the arc lengths, 


divided by their respective velocities, to arrive at an expression for time: 


asec, ‘ bsech, 


vy V2 


subject to the constraint that: 


atan(®,) + btan(®,) - C = 0 





Figure 4.3 Propagation of Light in Two Isotropic Homogeneous Media 
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Call the constraining function g(®,.®,), that is let: 


g(P,,®,) = atan(®,) + btan(®,)-C = 0 


By the method of Lagrange multipliers [Ref. 9], we form a new function h where: 


h = T- ag(®,,®,) 


dh ' 
Setting —— =—~ =——=0 yields the following three equations: 
, 0, OA 


seco tang, - hasec*, =0 
1 


7 sect tand, - Absec’, =0 
2 


atan ®, - btan ®,- C = 0 


which yield: 





asing, __ Aa 
v,cos*, cos’, 
bsind, Ab 





v,cos’, cos’, 
a tan ®,-btan ®,-C = 0 





We then have: 


Be aifcaag 
yy Vy 


which again yields Snell's Law: 


sind, sind, 


vy 16) 











3. Optically Anisotropic Media With Velocity a Function of Depth 

Now, with the help of the calculus of variations we consider the case 
where the velocity is not constant but rather a function of depth. In other words, 
consider the propagation of light in two optically anisotropic media where the 
velocity of the light particles is some function of y [Ref. 24] and [Ref. 26]. Again, 
using Fermat’s Principle, we arrive at the following equation which we want to 
minimize: 

T {2d - nVi-y? 
x, ds x Wy) 

This functional is of the form F(y,y’), therefore, it represents a special case of the 
Euler-Lagrange equation when the integrand has no x dependence, [Ref. 6] and [Ref. 


11]. Starting with the basic Euler-Lagrange equation: 
d oy! 
Pare oe =0 
We derive an expression which will simplify this second order differential equation 


into a first order differential equation: 
F-LP) =F -F y'-F,y" 
yg ee yyy EP yy® 


Then multiplying by y’: 


fad 
a) 


F,y'- Fy" -F,,y'y"= “IF -y'F,] 
Recalling that a first integral of a system of differential equations is a function which 
has constant value along the solution curves of the system, we can then say that the 
first integral of Euler’s equation is: 
F-y’Fy’ =C 
Using this equation, we now solve the differential equation: 


1 +y? . y'y’ agi 


MO) i+y? wo) 


From Figure 4.2, y’=tana, where >, ta Hence, 


sea tam’ 
seca(v(y)) seca(v{y)) 





or, 
cosa _G 
Vv 
therefore, 
cos( > -,) 
=C 
Vv 


from which we arrive at our familiar expression of Snell’s Law, mainly that: 


sin(>,) . sin($,) 


v 





V. 


1 2 


4. Optically Anisotropic Media With Velocity a Function of Direction 
Consider the propagation of light in two anisotropic media where the 

velocity of the light particles is some function of direction. That is v, = v,(q@) and 
Vv. = v,(8), where a and 8 are the angles formed by the incident ray and the 
horizontal, and the refracted ray and the horizontal, respectively, see Figure 4.4. We 


solve this problem using variational methods and show that the result is: 


1 od 
——[v,(B)sin(B)) = C 
vicp) 4B” 


which, for the simple case of constant velocities (homogeneous media) reduces easily 





1 d P 
—— Fo Vw rsin(a)] = 


v; (a) 


to our standard form of Snell’s Law: 


sing, _ sind, 
Y Mp 
Let y’ = tan(a), then a@ = tan’(y’). Applying Fermat's Principle. we wish to 


minimize: 




















Figure 4.4 Propagation of Light in Two Anisotropic Media 





pepe dx 


“1 vitan"'(y’) 


1 


SNe SRE THE fly). This is done 1o siniplify the arithmetic involved. Later, we 
v,(tan“'(y’)) 


1 


replace f(y’) by ——————_. 
v,(tan“"(y’4) 


Our problem is now to minimize: 
T= [*V1+y" fdr. 
x 
This functional is of the form F(y.y’), therefore we again make use of the special 
case of the Euler - Lagrange equation, [Ref. 6] and (Ref. 11]. 


Euler’s equation then is: 


F-y’Fy’ = Cor, 


foW1+y?-y'FO Wi +y?+ y" =o 0] =C 
l+y 





which further reduces to: 


LO yyyhfi+y? =C. 
1 +y? 
This by itself is a more general expression than that obtained by the method of 


Lagrange multipliers (Ref. 26]. Replacing our earlier definition of f(y’) and y’ we 


now have: 














—_—__—___——. - tana, 1 +tan 


v,(a)y 1 +tan 


cosa 


v,(a) 





vita a) 


1 +tan*a =C 


vi (a) 





-tanaseca 





or, 


=C. 


v,(a@) costa 








cosa sina ee 





vi (a) 


Now establishing a common denominator of v,*(a): 
v,(a)cose + sina vi(a) 
, ie 
Vv; (a) 
from which it can be seen that: 


1 








OF, 








We now desire to reach the same generalization of Snell's Law of 
Refraction using optimal control theory. 

We first include the derivation of the general problem, using optimal 
control. For a further explanation of this process, see Pontryagin (Ref. 11]. It ts 


ain twa points (x.v,) and (x,y,). see Figure 4.4, with a continuous, 


/ 
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piecewise smooth curve in such a way that the functional, Jy] = f "F(x,y,y Jdx reduces 
* 


to a minimum, where 


fay), if Gy)eM, 
Fanyy) = 
Gy) filzyy, if eM, 


define x° = [PF Gyy 46x! =x,x7 =y,y' =u 


From the maximum principle we have: 





dx°® dy, 
— =F J 9 ’ =0 
; (xy) es 
dc) _, ah _ OF 
: dt ox 
By at, _ OF 
dt : dt ay 


H = WF(xy.y’) + ¥, + hy’ 


To maximize H. we set 4 =0: 








—-= Yes +p, =0, where u = y’ 
Ou 


OF 
¥, FEF vox, 


Setting H = 0 we have: 
OF 
0 = WF y,u) +9, 5a 


YW, = —WoF (x, y,#) WSO 


From the jump conditions, ¥* = ¥ + yw grad g (x,y): 


ofr S00) ead neon 
dy’ dy’ 


Or, 
“h* = Ye¥, Zoran’ 
and. 
Yoo : 2 +N? 
oy oy 
or, 
oy ey 


where (N',N’) is a normal vector to the curve g(x,y) = 0 at the break point of the 
trajectory. Let Y’ denote the slope of the tangent to the curve of g(x,y) at that 


point: 





Be: ce 
pi tO ON 
Ff, = oY, +y/ / 
-f+—y y-2 yy’) Y 
a are arena 


from which: 
Or, 
f,+ + Seer (y*))=f+ ati (y)). 


Since the velocity of the light particles is a function of direction, the corresponding 


integral is: [PFO Vl+y" ae, where f(y’) = = 


Applying the results derived above, we have the following conditions: 


Vi OF a t35 
ot 

che aes 

ot oy 


Therefore, ¥, = C,, and ¥, = C,. Using this information, we now have: 


and 
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Ao dy’ 





We rewrite the right hand side of this equation as 


+y7- +y7 +4 y! 
AO Wi +y”? -lio Wi +y lice 


or, 
IO) -y fy V1 ry” 


l+y? 


= tan a and that: 











Recalling that y 
™ 
% = pani 
psi oy! de ety et 
fy) = Af, . a _ da _ sec’a oie 
gy) =— = 2 = = 
dy’ via)? via)? v,(a)? v,(a)? 
we then arrive at the expression 
dy, 
st # tanaseca 
v,(a) — v,(a)? 


which then becomes: 
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dv, 
v,(a)cosa +—sina 
da 
v,(a)? 


This is our familiar expression for this particular generalization of Snell's Law of 


Refraction: 


1d 
v,(a)? da 








[v,(o)sine] = 


Hplresinb. 


v,(B)? 


5. Optically Anisotropic Media Without the use of Lagrange Multipliers 
Consider again the same case where the velocity of the light particles 
depends upon the direction at each point. We wish to reach a generalization of 
Snell's Law without the use of Lagrange multipliers. A look at Figure 4.3 shows that 
we need to minimize the following expression: 
T's a Misses 
Ko) = -&(9,) 


where r = f(@.) is the polar equation of a smooth starshaped curve and r = g(®.) 





is similarly defined. A smooth starshaped curve has the property that every line 
segment connecting any point within the curve to the center of the curve remains 
completely within the curve. Another look at Figure 4.3 reveals that 

C = atan(®,) + btan(®,). 
Let h(®,.@,) = atan(®.) + btan(®,) - C. Recalling the definition of the total 


differential for a function, h, where h = h(®,,@.), we have: 
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oh oh 
dh = “1 do +“ ag 
b>, ' dh, 7% 


Then the following relationship holds: 


Now since we wish to determine the minimum of a function we simply set 


the derivative of the function equal to zero, that is: 


aT Xo. 
do, 


This, and the chain rule give us: 


[asece,tand (,) - asecd, f'(o,) ‘. 


| PO) 
bsecd,tand,e(,) > ~bsecd,g one 
1 1 a 4) 
g7(o,) 
which further reduces to: 


[aseco tand AY,) -asech,f(>,)| Z 


F,) 


oe isind,9(6,) ~ cOS,g8 (| 


8°(,) 
which yields: 


[Ao ,)sing, -f'(%,)cosd, | a [e(o,)sind, i 8'(,)cos,| 


F(o,) 2°(6,) 


This expression is yet another representation of Snell’s Law of Refraction for our 


given problem. 
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6. Optically Anisotropic Media in Two Dimensions 
Let us now consider the problem of the propagation of light in a two 
dimensional anisotropic media. That is, consider a light particle traveling from a 
point A to a point B. Let the velocities of the light particles be some function of 
two spatial variables, say y and z: v = v(y,z). For instance, suppose we arrange four 
glass containers, two on top of each other, and place different liquids in each 


container, resulting in varying particle velocities in each medium. See Figure 4.5. 





Figure 4.5 Propagation of Light in Two Dimensions 
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Applying Fermat's Principle, we arrive at our familiar expression using variational 


methods: 


T= Bds _ Byl+y?+2" 


Ad YA Wy,2) 
dt 


Our functional 


Jiy] = [Foy '22' dx 


has no x dependence. Therefore, starting with the two Euler - Lagrange equations 


for our situation: 
d 
ae - ay” =0 
and 
d 
Be at? =0. 


We again reduce these two second order differential equations as before 


and are left with: 


(F-y'F.). 


/ n nu, 4 
Py - Fyyy -Fyy¥Y ” de y 


Thus, in this case, Euler’s equation has the first integral: 
F-y’F, =C 


where C is a constant. Similarly for Z we arrive at: 


Applying these two Euler-Lagrange equations to our functional we see that: 
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vik Cones Bes 





[L+y? +27 , / C 
I 


y 
y = 
v(y,2) yl +y7 +27 v(y,z) 


and, 


Yl+y? +22 ak z! =C 
— ee 


~ 


V(y,2) V1 sy? +z? u(y,2) 
which yield: 


yl+y?+z” wy,2) 


and, 


L+y? 7 
vl+y? +z? wy) 


Cc 
Dividing the above two expressions and letting oe = C,we arrive at: 
2 


C,. 





Or, 


1 +2727 =C,( 1+ y’). 


If we let z’ = tan(a) and let y’ = tan(®) where 9, = 5 -a@ and >, = riba 


1 + tana) = C, (1 + tan’ (B)) 


or, 


sec(a) = C,sec’ (B) 








or letting C, = yc; 
sec (a) = C, sec (8) 


which in turn equates to: 


cosa = C, cosB where C, = x : 
4 


Remembering that 4, = 3 -a and >, = -fBwe now have the familiar form of 


Snell’s Law: 
sin(®,) = C, sin(®,). 


where C, represents the ratio of the velocities in the different media. 








V. APPLICATIONS 


A. INTRODUCTION 
In all of our previous derivations of Snell’s Law of Refraction, our functionals 
are of the form: 


Sy) = [flay +y? ae. 


This represents the integral of a function, f(x,y) with respect to the arc length s, 
where ds =1+y” dx. In this case, the Euler-Lagrange equation can be reduced to 


the following: 


HW 
f -fy’ ae Lae =0 
y L+y? 

where f = f(x,y). An explanation of this derivation is found in Gelfand and Fomin 
{Ref. 6:p. 19}. In our various developments of Snell’s Law, the function f(x,y), 
represents the reciprocal of the velocity of the light particles in question. Recalling 
from physics that: 

sind, v 

ascent Spe = My, [Ref.1] 

sind, v, 


we now introduce the term, index of refraction: n,,. Thus, in all of our earlier 


discussions we employ a functional which involves the index of refraction. In fact 


we have n(xy) = f(x,y). With this in mind, we now develop certain inverse 








problems to illustrate the importance of the index of refraction to various acoustical 


and optical problems. 


B. INVERSE PROBLEMS 
We pose the following problem: 
Given a certain host medium, find the index of refraction which will allow a particle 


of light to travel from a point A at (0,0) to a point B at (1,1) or to a point C at 


CG -1, 1) along the following curves: (1) a parabola (concave down), (2) a parabola 


(concave up), (3) a circle, (4) a cycloid, and (S) a straight line. See Figure S.1. 


1. Parabola 
We solve these different problems in order, starting with the parabolas. 
Since the parabola must have the following endpoints y(0) = 0, and y(1) = 1 we 


necessarily need the parabola: y = x’. Applying Euler's equation discussed above 


we arrive at: 








-n,y' = ny a 
l+y 
Or, 
-n,(2x) = 2n . 
1+4x 
A, 1 





ni x(1 +4x?) 


which, by partial fraction decomposition becomes: 


a2 








Tt /2-1 1 


(3) 


Figure 5.1 Particle Paths in an Isotropic Medium 
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This integrates to yield: 
In [n| = —In |x| Fin [1+4x2] +n IC}. 
Therefore, the desired index of refraction, is: 
Cyl +4? 
x 


The other parabola (concave up) also requires the same endpoint restrictions, 


n(x) = 


y(0) = 0 and y(1) = 1. Therefore, we force the equation of our parabola to be 


y=yx. Applying Euler’s equation, we are left with: 








s_ ty 
-nh —| 
n -1 
pe aise NE 
a gol 


or, 


oe eee 
“2¥x  yx(4x~1) 


from which it follows that: 





which, after integration becomes: 
in \n| = tin |4x+1) +In!C!. 
2 


So, the desired index of refraction for this parabola is the following: 
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n(x) = CV4x+1. 


2. Circle 
Now, we consider the circular path. Again, we desire that our curve pass 
through the points A and B. Therefore, our equation for the required circle is: 
x + (y-1 = 1 


which yields: 


y=yl-x? +1. 


From our earlier discussion, the Euler equation is: 





from which we get: 








Or, 





from which follows: 





Therefore, after integrating we arrive at: 


in |n| = -In |x| +In|C} =In|4} +n jc}. 
x 
Hence, the desired index of refraction for this medium, in order to travel! along a 


circle from point A to point B is: 


n(x) = ©. 
x 


3. Cycloid 
Recalling from differential equations, the parametric equations of a 


cycloid are of the form: 
k? , 
x= “7 (o-sind) 


k? 
Daal ZH ~cos) [Ref. 13]. 
Since we desire that our cycloid pass through the points A and C we necessarily 


reduce the above equations to: 


u 


x = @-sin® 
y = 1-cos®. 
By trial and error we determine that the function associated with a cycloid is only 


a function of y. As such, the corresponding Euler equation is: 








na = 
y L+y? 
Or, 
By 1 _2y" 
n 2 L+y? 


The chain rule and integration yield: 


[oredy = [Pw = [Dax 





AN A L+y? A 1+y? ax 
or, 
In'n| = 5 taj 1+¥4) +inic| 
so, 
ny) = Cyl +y”. 
a ; 
Now since 2 = @ =—Si9_ we have: 
dx & 1-cosh 
dd 
py J 
nyy.=C |e cm 
\ (1 -cosd) 
or. ArT ~ cosh). 
nly) = C | 21 cosh) © 
\ (.-cos@)? 


Therefore, the desired index of refraction for the cycloid is: 


y 


a 








4. Straight Line 
Our final problem involves the path of a straight line from point A to 


point B. That is, y = x. Using the Euler equation, we arrive at: ‘ 








n, = 
which integrates to yield the following: 
n(x) = C. 
Therefore, the index of refraction in the case of a straight line path is: 
n(x) = C. 
Computing the various indices of refraction does indeed give us the 
desired curve from a point A to a point B. For sake of brevity, we consider the first 
case, that of the parabola (concave down). That is, we show that: 


min Jp) = (Ves? 


2 
A - l+y" dk 


will yield the curve y(x) = x’. 


Applying the usual Euler-Lagrange equation, that is: 


d ‘ 
1 ee - P= 0, yields: 

















from which it follows that: 





y? cx? 
1 +y? 1 +4x? 
or, 
rn) o7x? 
1+(4-¢7)x? 
Setting C = 2 yields: 
y’ = 2x 


or. 
y(x) = x’ which is the equation for the parabola (concave down) passing from point 


A to point B in our problem. 


C. GENERAL EXPRESSIONS FOR THE INVERSE PROBLEM 


1. Function of X Only 
It is possible to derive a single, generic expression for a family of curves 
in terms of any function, f(x.y), where the various functions f(x,y) determine the type 


of curve generated. We now show the development of such an expression. We 


begin with our standard form of our functional: 





Ty} = f°feay V1 +y? ax. 


We first restrict our attention to those functions, f of the form f = f(x), that is, those 
functions that are only dependent upon x. This leaves us with: 
b 
Jy = f fey 1+y? de. 
Now, since the integrand does not depend on y, the corresponding Euler equation 


is: 





Fy =C 
or, 
foy' _¢ 
1 +y? 
which further reduces to: 
Pet HE 
[Ray -C? 


Now integrating the expression yields: 


cax 


yta= | ——— 
fix)? -c? 


This, then is the simple expression we have been looking for when the function f(x,y) 


is only dependent upon x. With this expression, we can now input any function and 
in turn generate a family of curves. As an example, we make use of an earlier 
problem, that of the path of a circle. In that development. we saw that the 


corresponding function for a circular path was: 











fix) = 


. 


#19 


For simplicity, we will let C = 1. That is, our problem involves fx) = and we 
x 


wish to show that this will generate a family of circles. 
Starting with the simple expression we have: 
cdx 


aD eee 


which becomes: 


y= |e 


[1 92 
\ 2 
Or, 


Letting u = 1-C°x’ we have: 


which becomes: 





1 
~ 2 7 eels 
(y-C, +x a 


which is the equation of a circle, centered at (0,C,) with radius 7 This 
Cc 


procedure can be followed for any function and will generate a family of curves 


associated with that particular function. 


2. Function of Y Only 
We now turn our attention to those functions which are dependent upon 
y alone. We derive an expression which will also produce a family of curves. As 
before, we start with the now familiar functional: 


Jiyl= {Ray i ry? ae. 


We consider only those functions dependent on y: 


Syl = f° +y? ae. 


Since the integrand has no explicit dependence on x, we employ the following Euler 


equation: 
F-y’F, =C 
So, we have: 
/ 
fy) [J+y?-y! Sy)y =C 
L+y? 
or, 
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R 
fy) -y41- LW =c 
L+y? 
which reduces to: 


fy) =Cyl+y” 


1 foY = Cc? 
y "Ge : 


After integration this becomes: 


or, 


xa =f _. 
ly -c? 


This then is the expression which produces any family of curves dependent only 
upon the function f(y). We now show that this expression, when given the proper 


function f(y), yields the family of curves for a cycloid. Recalling the earlier 


discussion, the corresponding function is f(y) = Cy? For simplicity we let C = 
2 


sl- 


This reduces the function to fly) = a Plugging this into the earlier derived 
y 


expression, we have: 





This leaves us with a Volterra integral equation of the first kind: 











x+@=C f ye Os : 
- ~C 
q 
For a further explanation and development of the Volterra integral equation, see 


Keener (Ref. 25:p. 101]. Using the Fundamental Theorem of Calculus and 


differentiating both sides with respect to y, we are left with: 








a&_C 
dy ; 
y-C? 
or, 
de C\ y-ct 
This reduces to: 
y? ry em a 


Or, 


| a aoe 
yy” +1) a 


This differential equation represents the equation of a family of cycloids [Ref. 27]. 


We include a solution to this differential equation for completeness. Let + =k? 
C 


and let y = k’sin’t. Then, we are left with the following expression: 
(1+y’)k’sin’*t = k° 


Or, 








This reduces to: 
y’ = cot(t). We consider only the positive square root in this case, since the slope 


of the graph is positive here. This then yields: 


dx = dy = tan(‘) dy. 


cot(t) 





However, dy = k’(2sin(t)cos(t))dt. Thus, 
dx = 2k’sin’tdt. 


We now let ® = 2t, apply the initial conditions, and make use of the double angle 


formula, sin*t = ,to arrive at: 


1 -cosd 
2 


k? 
dx = co (1 - cos) db 


are 
dy = 7 (sind) do . 


This yields: 


and, 


k? 
y= > (1 -cos) +C,. 


Using the fact that at @ = 0, x = 0 andy = 0 yields: 











2 
x= = ( -sind) 


and, 


k 
yey (1 -cosd). 
These are the equations of a family of cycloids. Again, this procedure can be 


followed for any function f(y) and will yield corresponding families of curves. 
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VI. CONCLUSIONS 


“" This thesis presents several derivations of Snell’s Law of Refraction and one 
of the principle of reflection. In all of these, we place different constraints upon the 
problem and produce the familiar laws governing reflection and refraction. The 
main result of this research is the expression found for Snell’s Law of Refraction in 
which the velocity of the light particles in the medium depends upon direction. We 
also use inverse problems to find the index of refraction needed to move along any 
desired family of curves. Finally, we derive a general expression that generates a 
family of curves given any index of refraction. 

.. This research has several applications in the area of acoustical and optical wave 
propagation. In the area of acoustics, the inverse problems enable us to predict the 
path followed by sound waves for a given medium. This would allow for military 
use in the area of prosecuting submarine contacts and in ocean bottom surveys. The 
inverse problems also have applications to the bending of light rays used in laser 
technology. 

Although this paper derives several generalizations of Snell's Law of Refraction 
there are still areas that require further research. For instance, consider the inverse 
problem in space where the index of refraction ts a function of three variables rather 


than two. Also, in all of the derivations in this paper the Hamilton-Jacobi equations 


p? 
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were not used. These equations, from the calculus of variations, could be used to 


derive a generalizaiion of Snell’s Law of Refraction. 
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